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ff-j \ Abstract The semi-classical back reaction to black hole evaporation (wherein 

the renormalized energy momentum tensor is taken as source of Einstein's 
equations) is analyzed in detail. It is proven that the mass of a Schwarzshild 
black hole decreases according to Hawking's law dM/dt = —C/M 2 where C 

— ! is a constant of order one and that the particles are emitted with a thermal 

q-i spectrum at temperature l/87rM(£). 
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On the basis of Hawking's |]J derivation of black hole radiation, wherein 
matter fields are quantized on the fixed background of a collapsing star, the 
mass of a black hole is expected to decrease according to the law 

dM C 

~df~~AP { > 

where C is a constant of order one which takes into account the number and 
spin of massless fields, the particles emitted at time t being distributed in a 
thermal spectrum with temperature 

T "W = sdm (2) 

Since Hawking's discovery much interest has focused on the gravitational 
back reaction in order to confirm or infirm (^) and (Q) and hopefully in the 
process learn something about quantum gravity. One of the simplest schemes 
in which to investigate the back reaction is the semi-classical theory wherein 
the renormalized energy momentum tensor is taken as source of Einstein's 
equations G^ v = 8Tr(T^ v ) ren , ie. the metric remains classical but follows self 
consistently the mean energy of the field. 

In this article we show that ([J) and (§) are indeed verified in the semi 
classical theory of black hole evaporation. However we mention at the outset 
that due to the exponential Doppler shift experienced by the field in its 
voyage from inside the star to X + the semi classical back reaction is not a 
consistent approximation to the true back reaction. This is not to say that 
the predictions of the semi classical theory are necessarily wrong, but the 
detailed mechanism by which vacuum fluctuations are converted into quanta 
will necessarily be totaly different in a fully quantum theory in which the 
metric as well as the matter fields fluctuate. 

This article is inspired by the work of Hajicek and Israel || and Bardeen 
[0] who showed that the semi classical back reaction is not incompatible 
with equation (|J) and by the work of Parentani and Piran |4| who showed 
numerically that (|I|) is verified in a model in which the renormalized energy 
momentum tensor takes a particularly simple form. 

We shall work in Bardeen's coordinates: 

ds 2 = -e 2i> {\ - 2m/r)dv 2 + le^dvdr + rW (3) 



In these coordinates Einstein's equations are 



4vrr 2 r; 



-Anr 2 T v v 



dm 

dv 

dm 

dr 

% = ^T rr (4) 

Note that ip is defined only up to the addition of an arbitrary function of v 
corresponding to a reparametrization of the v coordinate. 

We shall suppose that the right hand side of (^) is given by the renormal- 
ized energy momentum tensor (for a review see f|). The advantage of the 
coordinate system (||) is that the Hawking flux is encoded near the horizon 
in a negative energy flux T r v < 0. Were one to neglect the other components, 
the solution near the horizon would be a simple Vaidya type metric with only 
parameter m{y) = J v dv T vv . 

We shall proceed in three steps: first we shall suppose that the renormal- 
ized energy momentum tensor resembles the renormalized energy momentum 
tensor in the absence of back reaction. We shall then show that under this 
hypothesis the metric coefficients in (^) are slowly varying functions of r and 
v. Finally we shall solve adiabatically the Klein Gordon equation in this 
slowly varying metric and show that the renormalized energy momentum 
tensor indeed posseses the properties supposed at the outset thereby proving 
that the calculation is consistent. 

Our first task is to obtain estimates for T^ v both far and near the black 
hole. We begin with the former. We suppose that when r is equal to a 
few times 2m (say r = 0(6m)) there is only an outgoing flux T uu {r >> 
2m) = L^(w)/47rr 2 where Lh is the luminosity of the black hole. This is 
justified since in the absence of back reaction the other components of T^ u 
decrease as large inverse powers of r, for instance the trace anomaly decreases 
as m 2 /r 6 . We shall also suppose that Lh is small and varies slowly. Hence 
when r > 0(6m) an outgoing Vaidya metric is an exact solution of Einstein's 
equations 

ds 2 = _ (1 _ 2n^)) rfu 2 _ 2dudr + r 2 d Jl 

m(u) = J u du Lh(u) (5) 



The change of coordinates from (f|) to (||) is obtained by writing the equation 
for infalling radial null geodesies in the metric (||) as 

2dr 
Fdv = du+ — — (6) 

where F is an integration factor. Upon using (|6]) to change coordinates from 
the set (u, r) to (v,r) one finds that e^ = F and m(r,v) = m(u). Hence 
when r > 0(6m) the r.h.s. of (H) is given by 

A7rr 2 T r v (r > 0(6m)) = -e^L H 
-4nr 2 T v v (r > 0(6m)) = 2L H /(1 - 2m/r) 

4vrrT rr (r > 0(6m)) = AL H /r{l - 2m/rf (7) 

We now estimate T^ v near the horizon by assuming that the energy mo- 
mentum tensor measured by an inertial observer falling across the horizon 
is finite and of order L^. The jacobian relating v,r to (a),(/3) the or- 
thonormal frame of a free falling observer is e^ = (e~^(T — r) _1 ,r, 0,0), 
e*L = (e~^(r - r) _1 .r,0,0) where r = (1 - 2ffi + r 2 ) 1 / 2 and f is the 
rate of change of radius per unit proper time of the observer. If f ~ —1 
near r = 2m, the tetrad components are regular because T does not van- 
ish. Assuming T^^ = 0(L H /r 2 ) yields r 2 T r v {r ~ 2m) = 0(e^L H ), 
r 2 T v v (r ~ 2m) = O(Lh), rT rr (r ~ 2m) = OiLnjIm). The conservation 
of energy is 

(r 2 T;), r + r 2 T^ = (8) 

where r 2 T v vv = 0(Lh, v )- Integrating the conservation equation from r = 2m 
to r = 0(6m) yields T£ near the horizon in terms of its value where (0) is 
valid. Putting everything together, near the horizon we have 

— = -L H e^ + 0(mL H , v ) 

# = °^ 

#- = 0(L H /r) (9) 

o2m 

As anounced all metric coefficients vary slowly if Lh is small and varies 
slowly. Integrating the equation for ip yields e^ ~ r Lu for all r > 2m. Hence 



ip can safely be neglected up to distances r = 0(e 1 ^ LH ). From now on we 
suppose for simplicity of the algebra that ip — 0. 

Before proceeding we first discuss the location of the horizons in the 
geometry ([3]) with ip = 0. The apparent horizon is the locus where outgoing 
geodesies obey dr/dv = 0, therefore solution of 

r a (v) = 2m(v,r a (v)) (10) 

One may, following Bardeen, define the mass of the black hole at time v from 

M{v ) = r a {y)/2 whereupon (|9]) and (|I]) give 

^L = - Lh + o(L 2 h ) + 0(ML H , v ) (11) 

Hence equation ([[]) is recovered provided one proves Lh = C/M 2 and one 
takes v as time parameter. 

We also record here the equation for the event horizon th{v). This is 
the last light ray which reaches X + . It satisfies the equation of outgoing nul 
geodesies 

dr H _ l r H -2m(r H ,v) 

dv ~ 2 r H [ ' 

We obtain an assymptotic expansion for th{v) by requiring that it remains 
at a finite distance from r = 2M for all v , ie. neither diverges to r = oo nor 
falls into the singularity at r = 0. Setting rn{v) = 2M{v) + A(v ) we rewrite 
(HD in the form A = 2(2M+A)(2M tV + A >v ) + 2(m(v,M + A)-M). Solving 
recursively one obtains a series for A the first term of which is A = 8MM„. 
If one supposes L H = C/M 2 one finds that A = — 8C/M. 

In order to calculate the modes or (T^) we must first investigate the 
outgoing radial nul geodesies in the metric (|) with ip = 0. It is convenient 
to change variables to the set (v,x = r — r#(t>)) ie. x is the distance from 
the event horizon. In these coordinates the metric becomes (using (Ti2| ) ) 

ds 2 = 2m (v,r H + x)x dv2 _ ^^ _ r , dQ2 
ru\r H + x) 

When x « r H the equation for radial outgoing nul geodesies can be solved 
exactly to yield an exponential approach to the horizon of the form v—2 In x = 
f(u) where 

'"dv 2m ^ V)) (14) 

r 2 H {v) 



This motivates the following ansatz for the outgoing radial null geodesies 

2\nx + 5 = — T ^ + D (15) 



r H {v) 2m(u) 

with D a constant of integration. This should be compared with the solution 
in the absence of backreaction v — 2r — 4Mln(r — 2M) = u. The function of 
u on the r.h.s. of ([H|) has been writen as u/rh(u) for dimensional reasons. 
The quantity 5 is of order 0(Lh(Mx + x 2 )/M 2 ) for all x (this is shown by 



substitution of (|T5|) into the equation for radial nul geodesies and integrating 
the equation for 5 along the geodesies u = const). The function m(u) is 
determined by requiring that the variable u in equation (|15D be the same as 
in the Vaidya metric equation (|5|). The difference m(u) — m(u) is of order 
O(MLh)- This is found by using (|15"D to change coordinates from the set 
(f,r) to the set (u,r) at the radius r = 0(6M) where (^]) is valid. 

Equation ( |15| ) is sufficient to prove the first of our hypothesis, to wit that 
the flux emitted is 0(M~ 2 ). Indeed, neglecting the potential barrier and tak- 
ing only s- waves into account (as in the model of Parentani and Piran), the 
flux emitted is given by [|: T uu /Anr 2 = {l/\2ix)(dU / du) l l 2 (d/ du) 2 (dU / du)- 1 ' 2 
where U labels the outgoing geodesies as measured by an inertial observer 
inside the star (it is analoguous to the Kruskal coordinate U). The ja- 
cobian du/dU is calculated by remarking that at the surface of the star 
the derivative at fixed v dU / dr\ v=vstar = —2 is constant. Hence du/dU = 
— {\/2)du/dr\ v=vstar = (l/2)4m(u)/x\ v=vs tar = —4m(u)/U. In this expres- 
sion we have neglected transients which occur when the star has just col- 
lapsed (these transients vanish as e~ Au ^ im where Am is the time since the 
star began collapsing). Hence the flux at times An > 0(4m) is equal to 
(tt/12)T£(u) where T H (u) = l/8vrm(u) = l/8vrm(u)(l + 0(L H )) as an- 
nounced. (We mention that their is an intimate relation between the thermal 
flux T uu just calculated and the thermal distribution of the emitted photons. 
Indeed in the saddle point evaluation of Bogoljubov coefficients [[7|] the jaco- 
bian du/dU = —Arh{u) jU is the necessary ingredient to prove that the high 
energy photons are distributed Boltzmanly at temperature Th(u).) 

We now return to the second hypothesis and check that T^ u is regular near 
the horizon. To this end one writes the Klein Gordon equation in the coordi- 
nates (|T3| ) and notes that modes of the form <p\j tm = r~ 1 Y l m e~ tXv \x\ 2tX 6(±x) 
are solutions near the horizon |x| << M. The condition of positive frequency 
inside the star is described succintly (see j8| and ||) by asking that one may 
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continue x analytically in the upper half complex plane. This fixes the ratio 
of the 9(±x) pieces to be e _2?rA , hence the thermal character of the flux. 
When calculating Green's functions or derivatives thereof at the coincidence 
point one rescales A to A — > uj = A2m(t> ,r#)/r^ to obtain the same mode 
sum in the ultraviolet as in the absence of back reaction. Therefore the 
renormalized energy momentum tensor is regular near the horizon. This is 
sufficient for our purpose and proves that our scheme is consistent. 

One may however go beyond this qualitative discussion and make a quan- 
titative estimate of T^ u by using the following adiabatic expansion of the 
modes. One writes the solution of the Klein Gordon equation in the absence 
of back reaction as <^ m = r - l Y l m e- iu ' v \r-2M\ iiMu ' f(AMu, I, (r-2M)/2M). 
Then one makes the ansatz for the solution in the metric equation (|T^) 
VA,/,m — r~ 1 Y l m e~ lXv \x\ lX f(2X,l,x/r H (v))(l + x(x,v)). Calculation shows 
that x is proportional to Lfj and encodes effects such as a change of O(Ljj) 
in the potential barrier in the wave equation and a little frequency mixing 
in the potential barrier due to the time dependence of m(r, v). It gives rise 
to a modification of the energy momentum tensor which compared to the 
estimates in the absence of back reaction are ST^/T^ = O(Lh)- The valid- 
ity of this adiabatic expansion up to distances r = 0(6M) where the waves 
propagate as purely outgoing modes proves that all the emitted quanta, and 
not only the Boltzman ones, are distributed thermally with temperature 
Th{u) = l/8nm(u). 

In addition to the Hawking temperature Th being slightly different from 
l/87rm(w) their is also a slight athermicity due to the variation of the tem- 
perature over time. However none of these effects encode information about 
the quantum state of the star. In this respect the semi classical back reac- 
tion resembles the usual Hawking calculation in that one is confronted with 
a violation of unitarity for the outside observer fT0[|. 

We conclude with some remarks concerning the validity of the semi clas- 
sical back reaction used so far. In the preceding calculation we have made 
appeal to the geometry and the structure of the quantum field on scales ex- 



ponentially smaller than the Planck scale. This is manifestly absurd |TT| . We 
now investigate where in space time and at what scale will quantum back re- 
action effects completely modify the semi-classical back reaction. To this end 
we consider the fluctuations of T^ u and in particular the energy momentum 
correlated to the presence (or absence) of a specific Hawking photon at X + . 



We recall here the main properties derived in [[12] of the energy momentum 
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correlated to the presence of a particular particle (denoted io) at X + in a 
wave packet centered around uq, with energy ~ ujq = 0(M~ l ) : 1) it is given 
by the conserved tensor T£ u = (1/Anr 2 )f ln (u); 2) Jdu f io (u) = u ; 3) f l °(u) 
takes it values in a region of width u^ 1 around u ; 4) In addition to the 
piece of T*° just described in the outside region r > r H , their is a piece (the 
partner) on the other side of the horizon r < r#; 5) Inside the star T l ° takes 
the form (l/47rr 2 )/ l °(it) after reflecting at r = whereas before reflection it 
is the form 7*> = (1 / '4vrr 2 ) f° >(v) . 

To get a taste for back reaction effects we shall take T l ° as the source 
for Einstein's equations. This will give us a qualitative estimate of when 
back reaction effects become important. We first consider the interior of the 
star. Their one should express distances and energies in terms of the inertial 
coordinate U rather than u, hence reexpress T^° u as T$ u = T^° u (du/dU) 2 = 
T^(2m(w ) / 'x\ v=vstar ) 2 '■ The energy density (l/4:TTr 2 )T^ u reaches the Planck 
scale when r = 0(1) and x\ v=vstar ~ m(w )e~" 0//4m (" ) = 0(l)and the quan- 
tum field can then no longer be approximated as free in the center of the 
star. It is important to notice that the time necessary to reach this Planckian 
regime is Am = O(mlnm) which should be compared with the total evapora- 
tion time Am = 0(m 3 ). However this central region can perhaps be treated 
in a phenomenological way and could be irrelevant insomuch as Hawking 
radiation is considered. Hence we now consider the fluctuations of the geom- 
etry induced by the fluctuations of T^ v near the horizon. Then integration of 
equations (§) shows that the effect of the fluctuation T^ u = {1 / Anr 2 ) f l ° {u) 
induces a change in the mass parameter m — > m + 0(uq) in passing from 
u = u Q — 0{uJq 1 ) to u = u Q + OiujQ 1 ). However going backwards in time 
the photon Iq gets exponentially close to the horizon r# (since the light 
rays ( |i~5|) adhere exponentially to the horizon). It is very soon at a distance 
r — th = 0(M~ l ),ie. inside the domain where the horizon is shifting due to 
the presence of particle i itself. Hence the geometric concepts of apparent 
horizon, event horizon, outgoing light rays breaks down into a quantum soup 
(the quantum ergosphere of York |]14|) of width r — r^j = 0(M~ l ). In this 
region a free field approximation is certainly not valid. How the quantum 
theory of gravity contrives to build up (or maybe not to build up) Hawking 
radiation out of the soup is the subject of much conjecture. 
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